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A theory of intrinsic localized modes (ILMs) in anharmonic lattices is developed, 
which allows one to reduce the original nonlinear problem to a linear problem of 
small variations of the mode. This enables us to apply the Lifshitz method of the 
perturbed phonon dynamics for the calculations of ILMs. In order to check the the- 
ory, the ILMs in monatomic chain are considered. A comparison of the results with 
the corresponding molecular dynamics calculations shows an excellent agreement. 
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^ 1 Introduction 

■ It is already a well-known fact that localized vibrational excitations can exist 

in perfect anharmonic lattices [Tp|3P|5|6|7|8|9p0] . Such excitations are called 
intrinsic localized modes (ILMs) discrete breathers [S] or discrete solitons 
|1U] . The main approach, which has been used for the study of these vibrational 
excitations so far, is molecular dynamics (MD) simulations, which are based on 
the numerical integration of the classical equations of motion in a lattices with 
a finite number of the degrees of freedom. The approach is rather efficient in 
the case of simple ID lattices (see, e.g. [TT|12|13|14|15|16|17|18|19|2U|21|22] ). 
but requires very lengthy computations in the case of 2D and 3D lattices due 
to the rapid growth of the number of numerical operations with the increase of 
the number of the vibrational degrees of freedom (see [23|24II25] ). Therefore, it 
is of interest to develop other methods which would allow one, when studying 
ILMs, to reduce the amount of numerical computations. 
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Below we are developing a mean-field-type theory of ILMs which allows one to 
perform calculations for macroscopically large lattices of arbitrary dimensions. 
The theory is based on the consideration of small variations of the ILM am- 
plitude. The equations for these variations include characteristics of the ILM. 
This allows us to reduce the nonlinear problem of the ILM to a linear problem 
of the perturbed phonon dynamics; the perturbation comes from the ILM and 
it is determined self-consistently. To describe the effect of the perturbation we 
apply the Greens' function method of the perturbed local dynamics of phonons 
(the Lifshitz method), allowing us to calculate the ILM. As compared to the 
method which also uses the Greens functions of the perfect harmonic lat- 
tice for calculations of strongly localized ILMs, our method works generally, 
including the cases when ILMs have remarkable size. To test the theory, we 
consider the even and odd ILMs in a monatomic chain. In parallel, the molec- 
ular dynamics simulations of these ILMs are carried out. The comparison of 
the results of both calculations shows an excellent agreement. 



2 Self-consistency equations 

Let us start with the classical equations of the motion of atoms in a lattice 

MnUn = -Y. Kn!...n,_i«niWn2-.Mn,_i- (1) 

k nin2... 

Here Un is the Cartesian displacement of the atom n, M„ is its mass, V^'^\ V^^\ . 
are the harmonic and anharmonic springs; the subscripts n include both the 
site number and the number of the Cartesian component. We are consider- 
ing a localized solution of the equation ([1]) describing the stable ILM. The 
corresponding displacements of atoms are of the form 

Un{t) = C.n + AnCOSUlt + 0{uJl), 

where ui is the frequency of the ILM which lies outside the phonon spec- 
trum, \An\ is the amplitude of the ILM at the site n, ^„ is the shift of the 
equilibrium position of the atom n (i.e. the dc-component of the ILM; this 
component differs from zero due to odd anharmonicities), 0{ui) is the sum of 
the higher harmonics. The contributions of these harmonics are usually rather 
small [5|HI5|S|7f5|U] . Below they will be neglected. 

In our further consideration, we shall make use of the fact that M„(t + t)) is 
also a periodic solution of the equation ([1]). Taking r to be infinitesimal, we 
can find small variation of the ILM g„(t) = Un(t + t) — Unit) = —AnT smouit, 
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which also oscillate in time with the frequency ui and amplitudes ~ An (note 
the phase shift tt/2 of qn(t)); it satisfies the equation: 

MnQn = - J2i^nn' + <9Vanh/<9u„(9u„/ , (2) 
n' 

where Vanh is the anharmonic part of the potential energy. We consider the 
oc sin uit terms of this equation (i.e. we neglect higher-order harmonics of it). 
Taking into account that the oc sinojit term of the product cos"^" Uit sin Uit 
coincides with the time-independent term of the product 2 cos^" uit sin^ coit, 
we get 

Mnqn = -J2iVnn'+Vnn')qn', (3) 
n' 

where 

Vnn' = 2{sin^ Uit d'^Vs.nh/ dUndUn') . (4) 

Here the partial derivative is taken for Un = + cos uit, 

2tt/uji 


denotes the averaging over the period 27: /ui. 

The dc-shifts ^„ are not the independent parameters of the problem: they are 
determined by the ILM amplitudes. The corresponding relation reads [26] 

^n = Yl 9nn' (f^Vanh/ dUn') , (5) 
n' 

where w = G'J|JJ,(0)/VM„M„/, G'S|U,(0) is the static limit (w = 0) of the 
Green's function of the perfect lattice [27|28] 

Gnl'i.^) = ^njen'j/ (w^ - 
3 

Cnj OC e*'^^" is the polarization vector of the phonon j, kj is its wave vector. 
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Our further consideration is based on the observation that the equation ([3]) 
corresponds to the following harmonic potential energy 

= ^ H + Vnn')qnqn'- (6) 

nn' 



This potential energy, if the amplitude parameters An are chosen correctly, 
should lead to the appearance of the linear local mode, being an infinitesimal 
part of the ILM. The potential Vuf constitutes a mean field for the linear 
mode and, therefore for the ILM. The matrix f„„/ gives the required change 
of the elastic springs. 

The potential Vmf should be determined self-consistently. This can be done 
by applying the Lifshitz method [27] of the local dynamics in harmonic ap- 
proximation. In this method, the amplitude parameters satisfy the relations 

An/Ao = Gno{uJi)/Goo{iOi), (7) 



where G{u) is the matrix of the perturbed Green's function being equal 

G{uj) = {I- G'-'^^uj) (8) 



the frequency ui is given by the position of the pole of G{ijj) outside the 
phonon spectrum. Taking into account that the amplitude parameters An 
enter into the perturbation matrix v, one can conclude that equations ([7]) 
together with the pole condition (/ — G^'^\u)v)~'^ oo constitute the self- 
consistency equations. 

To calculate an ILM by means of the proposed method, one should choose 
the main amplitude (s) of the ILM and to fix a possible values of few other 
amplitudes. Using these amplitudes one can find the dc-shifts ^„ by solving 
the equation ([5]). Inserting obtained ^n into Eq. (jl]) one can calculate matrix v. 
After that one should calculate the frequency of the ILM and the relative am- 
plitudes, using equations ([8]) and ([7]). If the obtained An differ from the chosen 
initially, the calculations should be repeated with newly found amplitudes. 

The presented theory can be applied for different, including 3D, anharmonic 
lattices. 
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3 Test for monatomic chain 



As a test of the presented theory we consider the well-known case of ILMs in a 
monatomic chain, taking into account only the nearest-neighbor interactions 
and the positive quartic anharmonicity. The potential energy of vibrations in 
this case has the form 



V 



^ n 



K2 + ]^K^ul 



(9) 



where = — is the difference of the shifts of the particles n and 
n — 1. The top phonon frequency equals = 2y/K^. In our calculations we 
took K2 = 100 which corresponds to = 20. The Green's functions of the 
harmonic chain equal [2B] 

= 'l/c^v/^^^, (10) 

where p{u) = {u — y/ui'^ — 1)^ < 1 (in Eq. ( fTOl) the units = M = 1 are 
used) . 

The dc-shifts in this model equal zero, and the equation for the perturbation 
matrix v can be found strightforward: 

Vnn' = Sn,n'{ln+l + In) - 5n-l,n'7n " 5n+l,n'7n+l- (H) 



Here 



7n 



-K,At 



:i2) 



is the renormalization of the elastic spring between the atoms n and n — 1, 
An = An — An-i- FoT K4 > all 7„ are positive. Therefore the ILMs appear 
above the phonon spectrum; they have an odd or an even symmetry. 

First we study an even ILM centered at the atoms n = and n = —1. In 
this case, the amplitude parameters satisfy the condition An-i = —A^n- We 
consider only ILMs which have the largest amplitude on the atoms n = 
and n = —\ and take into account the contributions of 8 central atoms. 
In this case we include into consideration the renormalization of 7 central 
springs given by 70, 7±i, 7±2 and 7^3 (see Eq. (fT2|) ). These parameters can be 
determined together with the frequency and the amplitudes self-consistently 
using the equations ([8]) and ([7j). The latter we treat iteratively: we start with 
the amplitude ratios given in [S] and then find the frequency, the values of the 
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amplitudes and then their corrections from these equations and so on. The 
iteration procedure converges fast. The results of the calculations can be seen 
in Table 1 and in Fig. 1 (theory). 

Now we consider an odd ILM centered at the n = atom. Then the amplitude 
parameters of the ILMs satisfy the condition An = A^n- We restrict ourself 
with consideration only of ILMs having the largest amplitude on the central 
atom, and take into account the contribution of 7 central atoms. In this case 
we take into account the changes by the ILMs of 6 central springs given by 
70) 7±i! 7±2 and 73 (see Eq. f[T^ ). As in the case of even ILMs, the springs 
parameters are determined together with the frequency and the amplitudes 
self-consistently using the equations ([8]) and ([7]). The results of the calculations 
can be seen in Fig. 1 (theory) and in Table 1. 



4 Comparison with MD calculations 



In the case under consideration, one can easily find the ILMs numerically by 
integrating the equations of motion. In our calculations we used the forth- 
order Runge-Kutta algorithm. The results of calculations are compared with 
the theory. In Fig. 1 and in Table 1 one can see that the theoretical and the 
MD calculations of the ILMs frequencies and amplitudes are in very good 
agreement. 
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Fig. 1. Frequency dependences of the intrinsic localized modes in a monatomic 
chain with a hard quartic anharmonicity versus the dimensionless anharmonicity 
parameter K4Aq/4:K2] \Aq\ is the ILM amplitude on the central atom. 
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Table 1 

Frequencies and amplitudes of the even and odd ILMs in a monatomic chain with a 
hard quartic anharmonicity. For every K4 two values of the parameters mentioned 
are given: the theoretical value (up) and the result of the MD calculation (down). 



Even ILM 



K4 


K4 a2 




^0 




A2 


-^3 


50 


0.1250 
0.1224 


1.0963 
1.0945 


0.5000 
0.4947 


0.3211 
0.3208 


0.1537 
0.1550 


0.0662 
0.0680 


200 


0.5000 
0.4996 


1.3617 
1.3551 


0.5000 
0.4998 


0.1981 
0.1981 


0.0419 
0.0421 


0.0084 
0.0083 


800 


2.0000 

1.9999 


2.1289 

2.0995 


0.5000 

0.4999 


0.1218 

0.1219 


0.0084 

0.0085 


0.0005 

0.0005 



Odd ILM 



K4 


K4 a2 


1^1 


^0 




A2 


-^3 


50 


0.1250 
0.1210 


1.0862 
1.0839 


0.5000 
0.4920 


0.4034 
0.3989 


0.2266 
0.2288 


0.1054 
0.1117 


200 


0.5000 
0.4998 


1.3014 
1.2965 


0.5000 
0.4999 


0.3315 
0.3314 


0.0999 
0.1001 


0.0225 
0.0229 


800 


2.0000 
1.9920 


1.9228 
1.9040 


0.5000 
0.4990 


0.2858 
0.2726 


0.0386 
0.0358 


0.0031 
0.0029 



Note that the small discrepancies between the theoretical and the MD values 

of the ILM frequency Ui increases with Ui. The reason is that wc have neglected 
the higher harmonics: the contribution of these harmonics to an ILM is the 
larger the larger is the frequency. 



5 Conclusion 

To conclude, we have developed a theory of intrinsic localized modes in an- 
harmonic lattices, which allows one to make use of harmonic approximation 
results and to perform calculations for a macroscopically large lattice. We have 
derived the equations describing small variations of ILMs amplitudes, which 
allowed us to reduce the original nonlinear problem to the linear problem of 
phonon localization on the local effective potential. The latter is created by 
the ILM itself and it is determined self-consistently. This enabled us to apply 
the Greens' function method of the perturbed local dynamics of phonons (the 
Lifshitz method) for the calculations of ILMs. The theory works generally. 
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including the cases when ILMs have remarkable size. To test the theory we 
calculated even and odd ILMs in an anharmonic monatomic chain. We have 
also carried out the MD calculations of the ILMs in the chain. A comparison 
of the results, obtained by both methods (see Fig. 1 and Table 1) gives full 
confirmation of the proposed theory. Although concrete calculations here are 
performed only for the monatomic chain, the theory itself is general and can be 
applied for different lattices, including 3D crystals with arbitrary anharmonic 
potentials. 
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